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I.  INTRODUCTION 


It  is  now  widely  understood  that  poorly  designed  propelling  charges  can  lead  to  charge 
motion,  grain  fracture  (especially  for  cold  charges),  unprogrammed  burning  propellant  surface, 
high  chamber  pressures,  and  (occasionally)  breechblows.  Previous  studies  of  this  effect  using 
two-phase  flow  interior  ballistic  codes  have  been  performed  by  simulating  the  initial  action  of 
the  charge,  stopping  the  code  as  the  grains  are  shown  to  impact  the  rear  of  the  projectile, 
modifying  the  local  burning  surface  by  some  multiplicative  fector,  and  then  restarting  the 
calculation.  This  method  of  investigation  was  cumbersome,  but  it  yielded  an  improved  picture  of 
the  damage  that  could  be  wrought  by  an  ill-designed  charge,  e^)ecially  for  cold  ambient 
propellant  temperatures. 

In  the  present  study,  the  two-phase-flow  interior  ballistic  code  XNOVAKTC  (XKTC)  has 
been  modifi^  to  include  the  effects  of  propellant  grain  fracture.  An  increase  in  propellant 
surface  is  related  to  the  level  of  intergranular  stress  in  the  propellant  bed,  caused  either  by  local 
bed  compaction  associated  with  the  ignition  wave  or  by  grain  impact  against,  most  likely,  the 
proja:tile  base.  The  user  specifies  the  increase  in  propellant  surface  to  be  associated  with  the 
level  of  intergranular  stress.  An  increase  in  local  gas  production  follows  directly  from  the 
increase  in  burning  surface;  interphase  heat  transfer  and  drag  may  similarly  be  linked  to 
intergranular  stress.  This  improved  XKTC  is  used  to  illustrate  the  interior  ballistic  effects  of 
grain  fracture  for  two  charges  at  two  ambient  temperatures.  Pressure  waves  in  guns  have  long 
been  linked  to  localized  ignition  and  the  distribution  of  ullage  in  the  chamber;  propellant 
fracture  resulting  from  associated  increases  in  intergranular  stress  are  shown  in  this  study  to 
provide  a  link  to  higher  breech  pressures  as  well. 

II.  BACKGROUND 

Interior  ballisticians  have  long  studied  the  phenomenology  of  pressure  waves  in  guns.  The 
development  of  multiphase  flow  interior  ballistic  codes,  like  NOVA',  was  largely  motivated  by 
the  search  for  a  quantitative,  predictive,  and  diagnostic  tool  for  assessing  the  causes  and  cures 
for  pressure  waves.  These  studies  have  influenced  the  design  of  ignition  ^tems,  prop>ellant 
geometry,  and  even  charge  packaging  components.  By  and  large,  the  objective  has  b^n  to 
eliminate  or  minimize  pressure  waves,  for  overpressures  and  breechblows  have  long  been 
associated  with  gun  systems  which  exhibit  large  pressure  waves. 

The  problem  may  be  briefly  reviewed  in  reference  to  the  120-mm,  M829  APFSDS 
cartridge  depicted  in  Figure  1.  The  propellant  is  contained  in  a  combustible  cartridge  case 
which  is  affixed  to  the  projectile;  however,  the  propellant  often  does  not  fill  all  the  available 
volume  within  the  case,  particularly  after  some  settling  occurs  during  handling.  Further,  in  both 
tank  gun  and  artillery  charges,  the  ignition  system  provides  a  local  stimulus  which  leads  to  a 
nonuniform  ignition  of  the  main  propellant  charge.  During  the  ensuing  flamespreading,  the 
interplay  of  several  processes  determines  whether  pressure  waves  are  formed  and  their  influence 
on  the  overall  interior  ballistic  cycle. 

Once  ignition  occurs  at  some  point  within  the  charge,  there  is  competition  between  the 
local  gas  generation  rates,  determining  the  rate  of  local  gas  production  and  pressurization,  and 
macroscopic  charge  permeability,  determining  the  ease  with  which  this  local  accumulation  of 
gases  can  be  transported  to  the  rest  of  the  gun  chamber.  The  rate  of  gas  generation  is 
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controlled  by  the  intrinsic  burning  properties  of  the  prc^llant  as  well  as  by  its  geometry,  which 
influences  lx)th  surface  area  and  loading  density.  Macroscopic  charge  permeability  includes 
contributions  from  both  the  configuration  of  the  propellant  and  the  distribution  of  free  space,  or 
ullage,  within  the  chamber.  For  example,  the  natural  channels  through  a  bundle  of  stick 
propellant  offer  much  less  resistance  to  gas  flow  than  does  the  tortuous  path  through  a  bed  of 
granular  propellant.  Mobility  of  the  charge  itself  may  be  important,  as  early  dii^rsal  of  the 
propellant  bed  may  both  decrease  local  production  rates  and  increase  permeability. 
Unfortunately,  subsequent  impact  of  this  mobile  solid  phase  against  the  chamber  walls  or, 
especially,  against  the  projectile  base  can  lead  to  local  regions  of  high  loading  density  and 
propellant  fracture,  with  attendant  increases  in  burning  surface. 

Because  of  the  high  length-to-diameter  ratio  of  most  gun  chambers,  pressure  waves 
usually  manifest  themselves  in  the  longitudinal  mode.  In  high-loading-density  granular 
propellant  charges,  the  wave  finds  its  origin  in  the  pressure  gradient  that  forms  nearly 
coincidentally  with  the  convectively-driven  flame  front  traveling  from  an  ignition  source  in  the 
rear  of  the  charge  to  the  front  of  the  propellant  bed.  During  this  process,  the  pressure  gradient 
and  interphase  drag  forces  lead  to  the  development  of  intergranular  stresses  on  the  grains  and 
accelerate  them  to  impact  on  the  projectile  base.  A  reflected  compression  wave  in  the  gas 
phase  results,  whose  magnitude  may  be  increased  because  of  lower  porosity  due  to  bed 
compaction,  higher  gas  pressures,  and  especially  by  an  appreciable  increase  in  the  burning 
propellant  surface  area  increase.  Fortunately,  projectile  motion  works  to  expand  the  volume  for 
the  expanding  gas,  which  tends  to  attenuate  the  reflected  waves.  From  the  extensive  literature 
on  pressure  waves,  two  survey  references  are  cited 

Clearly,  long  bayonet  primers,  centercore  igniters,  and  other  distributed  ignition  systems 
all  represent  attempts  to  defuze  this  problem.  While  such  techniques  often  work,  their 
inadvertent  failure  (e.g.,  separation  of  a  primer  flashtube  from  its  headstock)  may  be 
catastrophic.  Further,  many  new  projectiles,  particularly  for  tank  guns,  protrude  far  back  into 
the  cartridge  case,  precluding  the  use  of  conventional,  long  bayonet  primers.  Coupled  with  the 
high  propellant  loading  densities  typifying  tank  ammunition,  some  level  of  longitudinal  pressure 
waves  is  to  be  expected. 


III.  NEWXNOVAKTC 

Like  all  versions  of  the  NOVA  code',  XNOVAKTC  (XKTC)  is  based  on  a  numerical 
solution  of  the  governing  equations  for  the  macroscopic,  quasi-one-dimensional  flow  defined  by 
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the  solid  propellant  and  its  products  of  combustion.  For  the  purposes  of  this  study,  the  principal 
extensions  of  NOVA  which  have  been  previously  been  incorporated  into  XKTC  are  the  ability  to 
model  the  intrusion  of  a  projectile  into  the  front  of  the  chamber  and  the  modeling  of  the 
ignition,  combustion,  and  compressibility  of  the  case.  The  intrusion  of  the  projectile  is  reflected 
in  the  cross  sectional  area  of  the  flow,  which  becomes  a  function  of  time  as  well  as  position. 

From  the  beginning,  the  NOVA  codes  have  recognized  the  role  of  propellant  bed  rheology 
as  significant  to  the  overall  interior  ballistic  process'.  While  individual  grains  are  still  treated  as 
being  incompressible,  compaction  of  an  aggregate  of  grains  is  allowed,  with  intergranular  stress, 
tr,  in  excess  of  ambient  pressure  taken  to  be  dependent  on  grain  porosity,  e,  and  the  loading 
history  of  the  bed.  The  constitutive  law  is  imbedded  into  the  formula  for  the  rate  of 
propagation  of  intergranular  disturbances,  a : 

J_ 
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where  is  the  density  of  the  solid  propellant,  and 
ga  is  a  constant. 

The  functional  dependence  of  intergranular  stress  on  fjorosity  is  shown  the  tne  Figure  2. 
Consider  first  the  upper  curve,  which  describes  the  stress  buildup  as  the  porosity  decreases  below 
the  settling  porosity,  e,,.  XKTC  assumes  that  the  speed  of  intergranular  disturbances  along  this 
curve  can  be  expressed  as 


«(0  = 


go  da 
Pr 


a(c)  = 


a,  €„ 


where  a,  is  the  value  of  a  when  the  bed  is  at  the  settling  porosity.  Thus  the  nominal  loading 
curve,  the  upper  curve  in  the  figure,  corresponding  to  monotonic  compaction  of  the  bed  from  e, 
to  a  smaller  value  of  porosity,  c,  is  given  by 


a  = 


I  2 
^0 


go 


Co 


When  the  bed  unloads,  or  reloads,  XKTC  assumes 


a(e)  =  a„ 

where  the  rate  of  propagation  of  intergranular  disturbances  in  this  unloading  bed,  Oj,  is  always 
larger  than  a,,  corresponding  the  much  higher  effective  modulus  during  the  unloading  process. 
This  results  from  the  effects  of  permanent  grain  deformation,  slippage,  and  grain  alignment 
during  loading.  Unloading  or  reloading  is  represented  on  the  figure  by  the  straight  lines  with 
large  negative  slope.  Values  of  a^  and  a,  are  deduced  from  laboratory  measurements  of  the 
stress  which  results  from  a  strain  imposed  upon  a  propellant  bed  . 

Consider  Figure  2  as  we  trace  the  process  of  bed  loading  and  unloading.  Scanning  the 
figure  from  right  to  left,  we  see  that  the  intergranular  stress,  a,  remains  zero  until  the  bed 
porosity,  e,  is  reduced  to  the  settling  porosity,  Then,  as  the  porosity  is  further  reduced,  the 
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Figure  2.  Intergranular  Stress  as  a  Function  of  Porosity 


stress  increases.  If  at  some  point  the  bed  is  unloaded  a  little  bit,  so  that  the  porosity  is 
somewhat  decreased,  the  stress  is  reduced  at  a  much  faster  rate  than  it  was  built  up,  according 
to  the  more  rapid  rate  of  change  of  stress  with  porosity  shown  in  the  figure.  If  the  stress  did  not 
reach  zero,  when  the  porosity  is  then  decreased  again,  the  stress  increases  rapidly,  back  to  the 
original  loading  curve.  The  loop  to  the  left  of  that  first  excursion  illustrates  the  effect  of 
unloading  the  bed  until  the  stress  is  reduced  to  zero  and  then  even  further  increasing  the 
prarosity.  In  that  case,  if  the  bed  is  again  loaded  by  decreasing  the  porosity,  the  stress  does  not 
follow  the  unloading  path,  but  instead  is  simply  assumed  to  increase  at  a  rapid  rate  of  change 
with  respect  to  porosity  to  the  original  loading  curve.  The  rates  at  which  loading  and  unloading 
take  place  are  characteristics  of  the  particular  propellant,  and  they  are  strong  functions  of  the 
propellant  temperature.  As  one  might  expect,  intergranular  stress  in  a  bed  of  cold  propellant 
builds  more  rapidly  than  it  does  in  propellant  at  ambient  temperatue. 

The  new  XKTC  has  certain  added  features  that  make  study  of  grain  fracture  easier. 
Without  manual  interference  during  the  calculation,  one  can  take  into  account  the  compaction 
of  the  propellant  bed  and  the  resulting  stress  on  the  propellant  grains  and  can  permit  the 
propellant  to  fracture,  increasing  the  gas  generation  rate,  interphase  drag,  and  heat  transfer  by 
an  amount  chosen  by  the  code  user.  The  extensions  to  XKTC  are  two-fold.  First,  constitutive 
modeling  has  been  extended  to  represent  the  influence  of  stress  induced  grain  fracture  on  the 
rate  of  gas  generation  and  the  interphase  drag  and  heat  transfer.  Second,  numerical  analysis  has 
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been  incorporated  to  treat  the  impact  induced  compaction  wave  traveling  backward  through  the 
bed  as  an  explicit  discontinuity. 

The  equations  of  motion  of  the  solid  propellant  reflect  the  stresses  created  the  contacts 
between  grains  as  well  as  those  of  gasdynamic  origin.  It  is  assumed  that  the  extent  of  grain 
fracture  at  any  position  in  the  propellant  bed  can  be  represented  as  a  function  of  the  maximum 
intergranular  stress  experienced  at  that  point  in  the  bed.  The  influence  of  fracture  is  assumed  to 
be  characterized  by  a  table  of  multipliers  whk:h  depend  on  the  maximum  stress.  Such  data  may 
be  acquired  by  placing  stressed  samples  of  propellant  in  a  closed  bomb  and  then  performing  an 
inverse  burning  rate  analysis  in  which  the  known  bum  rate  determines  the  total  burning  area. 
The  deoued  multiplier  is  then  determined  as  tlie  ratio  of  the  actual  area  to  that  predicted  by  the 
form  function  for  the  unstressed  propellant  At  each  time  step  during  the  calculation  the 
distribution  of  maximum  stress  is  scanned  at  each  mesh  point.  The  surface  area  of  the 
propellant  grains  is  determined  from  the  nominal  form  function  and  then  adjusted  in  accordance 
with  the  value  of  the  multiplier  corresponding  to  the  local  value  of  maxirnmn  stress.  The 
present  version  of  XKTC  does  not  include  a  provision  for  the  surface  area  multiplier  to  be  a 
function  of  the  fraction  of  the  grain  burned,  but  such  a  feature  could  be  added. 

The  code  permits  the  stress  history  to  affect  not  only  the  amount  of  burning  propellant 
surface  but  also  the  interphase  drag  and  heat  transfer,  since  these  are  also  sensitive  to  the 
surface  area.  In  some  cases,  however,  it  may  be  that  bum  rate  augmentation  occurs  as  a  result 
of  fractures  which  do  not  result  in  the  total  cleavage  of  the  grains.  Accordingly,  the  form 
function  multipliers  for  the  interphase  drag  and  heat  transfer  are  not  required  to  be  identical 
with  those  for  the  bum  rate  and  may  be  specified  independently. 

The  procedure  which  we  have  described  depends,  among  other  things,  on  a  reliable 
calculation  of  the  stress  field  throughout  the  propellant  bed.  Assuming  that  a  reasonable 
number  of  mesh  points  are  provided,  one  may  expect  the  calculation  of  the  stresses  created 
during  flamespreading  to  be  sufficiently  accurate.  However,  the  stresses  created  by  the  impact 
of  the  propellant  against  the  base  of  the  projectile  require  special  attention.  Since  the  boundary 
velocity  of  the  bed  may  drop  by  a  large  amount  at  the  instant  of  impact,  a  basic  assumption  of 
the  method  of  solution  is  violated  and  inaccuracies  may  be  expect^.  It  is  assumed  that  the 
state  variables  are  smooth  functions  of  position  and  time  so  that  their  derivatives  may  be  reliably 
approximated  by  finite  differences.  Clearly  this  will  not  be  true  at  the  front  of  the  bed  shortly 
after  impact.  Accordingly,  a  representation  of  the  impact- induced  shock  as  an  explicit 
discontinuity  has  been  incorporated  into  XKTC.  The  velocity  of  the  granular  shock  is 
determined  as  part  of  the  numerical  solution  until  the  shock  weakens  to  the  point  where  the 
jumps  in  the  porosity  and  grain  velocity  across  it  are  negligible.  The  methc^  of  analysis  is 
similar  to  that  previously  reported  in  connection  with  earlier  versions  of  the  NOVA  code.  It  was 
shown  that  while  the  errors  induced  by  a  failure  to  represent  the  granular  shock  explicitly  were 
negligible  as  far  as  predictions  of  maximum  pressure  and  muzzle  velocity  were  concerned, 
calculations  of  the  impact  induced  stress  levels  could  easily  be  in  error  by  100%.  Of  course,  the 
earlier  study  did  not  admit  the  possibility  of  grain  fracture.  If  it  had,  the  maximum  pressure 
would  not  have  proved  so  indifferent  to  the  granular  stress  field. 

IV.  BASELINE  DATABASE  FOR  THESE  CALCULATIONS 

We  chose  to  model  an  M829  cartridge  for  these  studies.  It  is  depicted  in  Figure  1.  Note 
that  the  case  is  combustible,  that  the  projectile  intrudes  significantly  into  the  chamber,  and  that 
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the  primer  is  thus  not  more  than  half  the  length  of  the  chamber.  The  basic  data  for  modeling 
this  round  comes  from  Robbins*,  who  performed  extensive  studies  to  achieve  a  reasonable 
model  for  the  combustible  case.  The  baseline  propellant  burning  rates  are  appropriate  for  JA2 
at  ambient  temperature.  For  all  of  the  calculations  reported  here,  analysis  of  the  explicit 
compaction  wave  in  the  solid  propellant  was  enabled  if  there  was  any  significant  impact  of 
propellant  on  projectile,  that  is,  if  the  relative  velocity  of  the  propellant  and  projectile  exceeded 
0.2^  m/s  (10  in./s).  For  the  instances  for  which  we  simulated  localized  ignition,  it  was  assumed 
that  the  primer  released  the  same  amount  of  energy  into  the  charge  in  the  same  ao^tion  time  as 
the  bayonet  primer,  but  that  all  of  the  energy  was  deposited  within  the  first  0.0254  m  (1  in.)  of 
the  ch^ge. 

We  obtained  grain  fracture  data  from  the  excellent  studies  of  high-rate  mechanical 
properties  of  propellants  by  Lieb’.  He  dox:uments  studies  of  M30  and  JA2  propellants,  so  his 
rqx»rt  includes  a  goxxl  description  of  the  behavior  of  the  JA2  contained  in  the  M829  cartridge. 
The  data  for  JA2  propellant  grains  that  have  been  damaged  by  firing  in  the  Gas  Gun  Impact 
Tester  (GGIT)  show  no  significant  fraorture  olamage  at  any  impaort  veloxiity  within  the  range  of 
the  tests  if  the  temperature  of  the  grains  was  higher  than  -10°C.  For  propellant  temperatures  of 
-20°C  and  below,  however,  damage  did  oxxur,  and  it  was  quite  significant  for  impact  velox:ities 
of  90  and  115  m/s.  These  tests  show  that  impact  of  cold  grains  against  a  solid  target  at 
velocities  of  a  few  tens  of  meters  per  second  can  lead  to  surface  area  increases  of  up  to  a  factor 
of  twelve  (12).  It  is  interesting  to  note  that  the  surface  area  ratio  (surface  area  after  fracture 
divided  by  original  surface  area)  is  a  strong  function  of  the  fraction  of  the  grain  burned;  the 
highest  multiplicative  factors  are  associated  with  the  initial  burning  of  the  grain.  The  data 
clearly  show  that  if  the  grain  is  cold  enough,  impact  velocities  of  60,  90,  and  115  m/s  result  in 
sufficient  stress  of  the  grain  to  cause  significant  fracturing.  One  is  thus  able  to  associate  the 
stress  caused  by  ignition  and  especially  impact  with  the  stress  level  necessary  for  fracture.  We 
have  not  used  Lieb’s  data  directly,  as  aggregate  data  could  be  used  somewhat  more  directly  in 
the  present  version  of  XKTC.  We  have,  however,  used  the  data  as  a  guide  as  we  performed  the 
parametric  studies,  as  will  be  discussed  later. 

Two  other  propellant  mechanical  properties  should  be  considered.  First,  Lieb’  has  shown 
that  JA2  flows  plastically,  and  this  occurs  at  higher  values  of  stress  as  the  temperature  decreases. 
XKTC  includes  the  effects  of  propellant  deformation  only  as  it  is  involved  in  the  reduction  of 
the  porosity  as  the  propellant  bed  is  compressed.  Second,  Costantino  and  Omellas’  have 
shown  that  the  higher  the  gas  pressure  surrounding  JA2,  the  tougher  the  propellant  is.  We  have 
not  included  this  effect.  The  observation  of  Costantino  and  Omellas  provides  another  reason  to 
pressurize  gun  propellant  charges  evenly,  to  preclude  the  possibility  of  an  intergranular  stress 
wave  moving  ahead  of  the  pressurization  wave  into  weaker  propellant. 

The  resulting  baseline  data  base,  and  the  XKTC  echo,  is  included  as  Appendix  A.  Note 
that  the  front  of  the  combustible  case  has  been  given  a  tapered  front  end,  to  minimize  problems 
that  can  result  with  sudden  area  changes  in  a  one-dimensional  code  like  XKTC. 

V.  CALCULATIONS 

The  first  simulation  was  for  the  M829  cartridge,  with  its  standard  bayonet  printer,  at 
ambient  temperature.  Figure  3  shows  the  simulated  pressures  at  the  breech,  at  the  front  of  the 
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Figure  3.  Pressures  and  Pressure  Difference  for  Baseline  Case 
Ambient  Temperature,  Bayonet  Ignition,  No  Front  Ullage,  No  Propellant  Fracture 


chamber,  and  their  difference.  The  maximum  breech  pressure  for  this  case  is  508.7  MPa,  and 
19.2  MPa  is  the  first  negative  pressure  difference.  As  discussed  before,  XKTC  keeps  track  of 
the  maximum  stress  experienced  at  each  location  in  the  propellant  bed.  Figure  4  shows  the 
maximum  stress  for  this  simulation,  plotted  just  after  all  the  stresses  in  the  bed  had  relaxed. 
While  there  was  no  impact  of  propellant  on  the  base  of  the  projectUe  in  this  calculation,  the 
ignition  wave  is  seen  to  have  pressed  the  propellant  up  into  the  diminishing  area  behind  the 
projectile,  resulting  in  markedly  increasing  maximum  stress  at  the  front  of  the  bed.  The 
maximum  stress  reached  at  the  front  of  the  bed  is  26. 1  MPa.  Note  the  small  stress  at  the  rear 
of  the  bed,  which  should  be  expected,  as  this  bayonet  primer  causes  the  bed  to  be  driven  in  both 
directions.  Figure  5  shows  the  progression  of  intergranular  stress  through  the  bed  as  a  function 
of  time.  Each  of  the  stress  curves  has  been  offset  to  permit  one  to  see  the  stress  wave.  Time 
advances  from  front  to  rear  of  this  figure  at  0.1  ms  per  curve.  Note  the  stress  buildup  against 
the  projectile  and  then  the  rapid  decay  of  the  stress  as  the  bed  unloads. 

Figure  6  shows  the  pressures  that  result  if  one  permits  the  stresses  induced  in  the  first 
simulation  to  fracture  the  propellant.  The  stress  versus  propellant  surface  area  multipliers  that 
were  used  for  this  second  simulation  are  shown  in  Table  1. 
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Figure  4.  Maximum  Intergranular  Stress  in  the  Propellant  Bed 
Ambient  Temperature,  Bayonet  Ignition,  No  Front  Ullage,  No  Propellant  Fracture 
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Figure  5.  Intergranular  Stress  Versus  Position  in  the  Bed  as  a  Function  of  Time 
Ambient  Temperature,  Bayonet  Ignition,  No  Front  Ullage,  No  Propellant  Fracture 
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Figure  6.  Pressures  and  Pressure  Difference  for  Second  Simulation 
Ambient  Temperature,  Bayonet  Ignition,  No  Front  Ullage,  With  Propellant  Fracture 


Table  1.  Propellant  Surface  Area  Multipliers 
Propellant  at  Ambient  Temperature 


Stress  (MPa) 

0.0 

10.0 

25.0 

50.0 

Multiplier 

1.0 

1.0 

1.5 

4.0 

Between  stress  values,  the  multiplier  varies  linearly.  Above  the  maximum  stress,  the 
multiplier  remains  constant  at  the  maximum  value.  As  one  can  observe  from  Figure  5,  stress  of 
this  magnitude  occurs  only  at  the  very  front  of  the  bed,  due  to  compaction.  We  see  that  the  use 
of  these  data  leads  to  sufficient  grain  fracture  to  cause  the  maximum  breech  pressure  to  increase 
to  521.7  MPa,  and  the  first  negative  pressure  difference  increases  to  25.45  MPa,  as  shown  in 
Figure  6.  Figure  7  plots  the  maximum  intergranular  stress  in  the  bed  for  the  second  simulation, 
and  Figure  8  shows  the  time  progression  of  the  stress  wave  in  the  bed.  The  maximum 
intergranular  stress  produced  in  this  simulation  is  25.4  MPa. 

The  third  simulation  is  for  a  cold  charge.  The  propellant  burning  rates  have  been  adjusted 
to  be  more  appropriate  for  cold  temperatures*.  The  speed  of  the  Intergranular  disturbances  in  a 
settled  bed  has  b^n  increased  from  254  m/s  (10000  in./s)  at  ambient  temperature  to  381  m/s 
(15000  in./s).  The  database  for  this  third  simulation  is  included  as  Appendix  B.  The  pressures 
and  pressure  difference  for  this  simulation  are  shown  in  Figure  9.  As  expected,  due  to  the  lower 
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Figure  8.  Intergranular  Stress  Versus  Position  in  the  Bed  as  a  Function  of  Time 
Ambient  Temperature,  Bayonet  Ignition,  No  Front  Ullage,  With  Propellant  Fracture 
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Figure  9.  Pressures  and  Pressure  Difference  for  Cold  Charge 
Cold  Temperature,  Bayonet  Ignition,  No  Front  Ullage,  No  Propellant  Fracture 

propellant  burning  rates,  the  pressures  are  all  lower  than  those  observed  for  the  ambient 
calculations;  the  maximum  breech  pressure  is  380.2  MPa,  and  the  first  negative  pressure 
difference  Ls  14.0  MPa.  On  the  other  hand,  due  to  the  faster  wave  speed,  the  intergranular 
stress  is  increased  for  this  cold  case,  as  shown  in  Figure  10.  The  maximum  intergranular  stress 
produced  in  this  simulation  is  43.2  MPa.  The  time  progression  of  the  stress  wave  for  this 
situation  is  shown  in  Figure  11. 

The  fourth  simulation  is  for  a  cold  charge,  but  the  propellant  has  been  permitted  to 
fracture  according  to  the  multipliers  shown  in  Table  2.  As  lieb  found’,  there  is  more  fracture 
at  colder  propellant  temperatures. 

Table  2.  Propellant  Surface  Area  Multipliers 
Cold  Propellant 


Stress  (MPa) 

0.0 

10.0 

25.0 

0 

50.0 

Multiplier 

1.0 

1.0 

2.0 

10.0 

The  pressures  and  pressure  difference  for  a  cold  charge  with  propellant  fracture  are  shown 
in  Figure  12.  The  peak  pressure  is  up  to  427.1  MPa,  and  the  first  negative  pressure  difference  is 
now  37.5.  The  maximum  intergranular  stress  in  the  cold  bed  is  shown  in  Figure  13,  and  the  time 
progression  of  intergranular  stress  is  shown  in  Figure  14. 
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Figure  10.  Maximum  Intergranular  Stress  in  the  Propellant  Bed 
bid  Temperature,  Bayonet  Ignition,  No  Front  Ullage,  No  Propellant  Fracture 


Figure  11.  Intergranular  Stress  Versus  Position  in  the  Bed  as  a  Function  of  Time 
Cold  Temperature,  Bayonet  Ignition,  No  Front  Ullage,  No  Propellant  Fracture 
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Figure  12.  Pressures  and  Pressure  Difference  for  Cold  Charge 
Cold  Temperature,  Bayonet  Ignition,  No  Front  Ullage,  With  Propellant  Fracture 


Figure  13.  Maximum  Intergranular  Stress  in  the  Propellant  Bed 
Cold  Temperature,  Bayonet  Ignition,  No  Front  Ullage,  With  Propellant  Fracture 
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Figure  14.  Intergranular  Stress  Versus  Position  in  the  Bed  as  a  Function  of  Time 
Cold  Temperature,  Bayonet  Ignition,  No  Front  Ullage,  With  Propellant  Fracture 


The  maximum  intergranular  stress  produced  in  this  fourth  simulation  is  39.1  MPa,  which  is 
actually  smaller  than  that  produced  when  there  is  no  propellant  fracture.  To  repeat,  while  the 
difference  is  small,  the  maximum  stress  for  the  case  involving  grain  fracture  is  actually  smaller 
than  the  comparable  case  when  fracture  is  disabled.  Very  close  examination  of  the  details  of  the 
flow  shows  that  this  counterintuitive  result  can  be  explained.  For  the  first  1.3  ms  of  the 
simulation,  conditions  in  propellant  bed  A  (no  fracture)  and  propellant  bed  B  (fracture  enabled) 
are  essentially  identical.  At  about  that  time,  the  level  of  the  stress  begins  to  cause  grains  to 
fracture  in  bed  B.  The  grains  that  fracture  do  so  at  the  front  of  the  bed,  so  the  gas  pressure 
builds  faster  in  the  front  of  bed  B  than  that  in  bed  A.  This  localized  pressure  buildup  in  the 
front  of  bed  B  causes  the  flow  of  gas  toward  the  projectile  to  be  reduced,  which  in  turn  causes 
the  propellant  grain  velocities  at  the  front  of  bed  B  to  be  reduced.  For  a  few  0.01  ms  time 
steps,  in  fact,  the  net  gas  flow  at  a  few  stations  in  bed  B  is  away  from  the  projectile,  and  the 
grain  velocities  (toward  the  projectile)  there  and  farther  back  in  the  charge  are  noticeably 
smaller  than  those  in  bed  A.  As  a  result,  the  dynamic  conditions  in  bed  B  cause  the  peak  stress 
that  is  reached  to  be  somewhat  less  than  that  reached  in  bed  A,  in  which  no  localized 
gasdynamics  worked  against  the  motion,  packing,  and  stressing  of  the  grains. 

This  effect,  while  small,  results  from  very  tightly  coupled  physics.  Stress  leads  to  fracture, 
which  leads  to  gas  generation,  which  affects  gas  motion,  which  affects  grain  motion,  which 
inhibits  stress  build-up. 
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Serious  problems  could  exist  if  localized  ignition  were  to  take  place,  such  as  could  happen 
if  the  bayonet  primer  fractured  at  its  base,  so  that  all  of  its  energy  were  released  very  close  to 
the  breech  of  the  chamber.  The  situation  would  be  even  worse  if  some  propellant  settling  had 
taken  place,  creating  significant  front  ullage.  Such  a  situation  might  be  obtained  if  tank  rounds 
were  stored  vertically,  projectile  up,  if  the  propellant  settled,  and  if  it  did  not  immediately 
recover  to  its  original  cartridge-filling  configuration  on  loading.  TTiis  combination  of  unfortunate 
circumstances  is  simulated,  with  grain  fracture  permitted,  in  the  fifth  simulation.  The  database 
for  this  simulation  is  included  as  Appendix  C.  The  pressures  and  pressure  difference  are  shown 
in  Figure  15.  The  maximum  breech  pressure  is  "onl^  579.1  MPa,  but  the  first  negative  pressure 
difference  is  now  197.3  MPa.  Figure  16  shows  the  maximum  stress  for  this  simulation,  plotted 
just  after  all  the  stresses  in  the  bed  have  relaxed.  Not  only  is  there  profound  stress  on  the  front 
of  the  bed,  due  to  compression  and  impact,  there  has  been  significant  stress  even  on  the  center 
of  the  bed,  due  to  the  wave  from  the  localized  ignition  source.  The  relative  velocity  of  the 
propellant  and  the  projectile  at  the  moment  of  impact  (the  "impact  velocity")  was  about  75  m/s. 

The  intergranular  stress  versus  bed  position  as  a  function  of  time  is  shown  in  Figure  17. 
Note  the  greatly- increased  stress  in  the  center  of  the  bed  that  results  from  the  localized  bed 
ignition.  One  can  observe  that  the  stress  wave  grows  in  strength  as  it  moves  toward  the 
projectile,  but  that  it  collapses  to  zero  when  the  ullage  is  reached,  and  the  bed  unloads.  The 
maximum  intergranular  stress  produced  in  this  simulation  is  53.2  MPa. 


Time  (ms) 

Figure  15.  Pressures  and  Pressure  Difference  for  Cold  Charge  with  Problems 
Cold  Temperature,  Local  Igiition,  Front  Ullage,  Propellant  Fracture 
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Figure  16.  Maximum  Intergranular  Stress  in  the  Propellant  Bed 
Cold  Temperature,  Local  Ignition,  Front  Ullage,  Propellant  Fracture 


Figure  17.  Intergranular  Stress  Versus  Position  in  the  Bed  as  a  Function  of  Time 
Cold  Temperature,  Local  Ignition,  Front  Ullage,  Propellant  Fracture 
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Figure  18  shows  the  details  of  the  stress  history  around  the  time  of  maximum  stress.  Here 
the  curves  are  only  0.02  ms  apart.  One  sees  that,  after  the  stress  had  collapsed  when  the  peak 
reached  the  ullage,  it  later  builds  again,  well  away  from  the  projectile,  even  before  the  moving 
bed  impacts  the  rear  of  the  projectile.  Very  close  examination  of  the  details  of  the  calculation 
reveals  the  details  of  this  unexpected  finding: 

•  Igniter  functions.  Bed  ignites  from  the  rear.  Rames  reach  front  of  bed. 

•  Forward  motion  of  gas  causes  forward  motion  of  grains;  stress  wave  is  observed 

moving  forward. 

•  Gas  pressure  builds.  Projectile  starts  to  move. 

•  Wave  of  gas  arrives  at  projectile,  flows  across  ullage  in  front  of  grains. 

•  Gas  rebounds  off  projectile,  so  gas  now  opposes  grains  rushing  toward  projectile. 

•  Gas  pushes  back  on  grains,  builds  porosity  minimum  and  stress  maximum  well  away 

from  the  projectile  surface. 

•  Before  grain  impact,  front  grains  are  slowing  and  projectile  is  accelerating. 

•  At  the  time  of  impact,  there  is  a  solid  stress  peak  well  away  from  projectile. 

•  Impact  causes  sharp  local  drop  in  porosity,  sharp  rise  in  stress. 

•  Compaction  wave  begins  to  fiU  in  the  stress  curve  behind  impact  peak. 

•  Gas  velocity  strongly  opposes  grain  flow,  slowing  grains. 

•  Compaction  wave  dies  out  --  bed  unloads. 

•  Stress  rapidly  declines  to  zero. 


Figure  18.  Intergranular  Stress  Versus  Position  in  the  Bed  as  a  Function  of  Time 
Details  About  the  Maximum,  Curves  Each  0.02  Millisecond 
Cold  Temperature,  Local  Ignition,  Front  Ullage,  Propellant  Fracture 
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VI.  DISCUSSION 


We  found  that  the  interior  ballistic  implications  of  grain  fracture  for  conditions  which 
involve  cold  propellant,  localized  ignition,  and  front  ullage,  singly  or  in  combination,  are 
profound.  These  conditions  are  sufficient  to  lead  to  large  pressure  waves  and  greatly  enhanced 
peak  pressures.  The  most  important  and  damaging  propellant  stress  results  from  propellant 
impact  on  the  projectile,  but  the  stress  on  the  prop>ellant  bed  from  the  rapid  ignition  wave  can 
feed  a  potentially  dangerous  pressure  wave. 

We  used  Lieb’s  data’  in  the  following  manner.  The  use  of  the  gas  gun  impact  tester  and 
the  miniature  closed  bomb  produced  plots  of  the  surface  area  ratio  (the  ratio  of  the  surface  of 
fractured  grains  to  the  surface  of  an  unfractured  grain)  versus  the  fraction  of  the  grain  burned, 
as  a  function  of  the  impact  velocity,  for  several  temperatures.  We  noted  that  for  an  impact 
velocity  of  a  few  tens  of  meters  per  second,  especially  for  cold  ambient  temperatures,  the  area 
ratios  determined  could  range  from  one  to  twelve,  with  higher  impact  velocities  and  colder 
temperatures  resulting  in  higher  ratios.  In  test  calculations,  we  observed  that  for  local  ignition, 
front  ullage,  and  no  grain  fracture,  the  predicted  impact  velocity  was  several  tens  of  meters  per 
second  (specffically,  80  m/s),  so  we  chose  levels  of  stress  and  values  of  surface  multipliers  that 
gave  sig^iificant  grain  fracture  under  those  circumstances.  We  are  aware  that  lieb  plans  to 
acquire  fracture  data  for  a  propellant  aggregate  We  support  the  need  for  such  data,  for  they 
will  describe  very  directly  the  effects  of  compaction  on  the  prop)ellant  in  a  gun  chamber.  We 
look  forward  to  being  able  to  use  the  new  data  in  subsequent  studies. 

VII.  CONCLUSIONS 

XNOVAKTC  has  been  modified  to  permit  studying  the  interior  ballistic  effects  of  grain 
fracture.  This  new  capability  applies  not  only  to  grain  fracture  caused  by  the  impact  of 
propellant  grains  on  axid  boundaries,  such  as  the  projectile,  but  also  to  fracture  produced  the 
ignition  wave.  Further,  the  intergranular  stress  in  the  propellant  bed  caused  by  the  interaction 
of  that  wave  with  the  reduction  in  cross  sectional  area  due  to  projectile  intrusion  leads  to 
increased  fracture.  All  of  these  capabilities  are  folly  incorporated  into  the  code,  so  the  many 
interactions  of  physics  and  chemistry  involved  in  the  interior  ballistic  cycle  may  be  studied 
together. 

The  studies  included  in  this  report  support  earlier  findings  that  propellant  grain  fracture 
may  lead  to  increased  pressure  waves  and  increased  breech  pressures  for  cold  rounds,  especially 
if  there  is  localized  ignition  of  the  propellant  bed  and  there  is  significant  axial  ullage.  The  flow 
characteristics  of  stick  propellants,  which  minimize  their  motion  in  the  gun  chamber,  certainly 
recommend  their  use  in  an  attempt  to  minimize  circumstances  which  lead  to  prain  fracture. 
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APPENDIX  A 
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XNGVMCTC  VERSION  NUMBER  1.009 

CONTROL  DATA 

LOGICAL  VARIABLES: 

PRINT  T  PLOT  F  DISK  IRITE  F  DISK  READ  F 
I.M.  TABLE  T  FLAME  TABLE  T  PRESSURE  TA8LE(S>  T 
EROSIVE  EFFECT  0  IMLL  TEMPERATURE  CALCULATION  0 
BED  PRECOMPRESSa)  0 
HEAT  LOSS  CALCULATION  1 

BORE  RESISTANCE  FUNCTION  1 

SOLID  TRAVELING  CHARGE  OPTItM  (0=M0.1::VES)  0 

EXPLICIT  COMPACTION  UAVE(0=NO;1=YES)  1 

CONSERVATIVE  SCHEME  TO  IHTEGRATE  SOLID-PHASE  CONTINUITY  EQUATION  (O-NO.OLO;  IsTES.NEU)  0 
KINETICS  NODE  (0=HaNE;1=GAS-PHASE  0MLY;240TH  PHASES)  0 
TANK  GUN  OPTION  (0=N0.1=YES)  1 

INPUT  ECHO  OPTION  0 

FORWARD  BOUNDARY  CONDITION  (0=CLOSED;1^XiEH;2=ROCKET)  0 


LIQUID  TRAVELING  CHARGE  0PTiaN(0=N0,1=TES)  0 

GRAIN  FRACTURE  QPTI0N(0=NO;1=YES)  1 

GRAIN  FRACTURE  DATA  BASED  ON  INTRINSIC  AVERAGE  STRESS 

(0=N0;1=YES)  1 

INTEGRATION  PARAMETERS 

NUMBER  OF  STATIONS  AT  WHICH  DATA  ARE  STORED  99 

NUMBER  OF  STEPS  BEFORE  LOGOUT  3500 

TINE  STEP  FOR  DISK  START  0 

NUMBER  OF  STEPS  FOR  TERMINATION  3500 

TIIC  INTERVAL  BEFORE  LOGOUTISEC)  0.1000E-OL 

TINE  FOR  TERMINATION  (SEC)  0.4500E-02 

PROJECTILE  TRAVEL  FOR  TERMINATION  (INS)  186.90 

MAXIMUM  TIME  STEP  (SEC)  0.1000E-03 

STABILITY  SAFETY  FACTOR  2.00 

SOURCE  STABILITY  FACTOR  0.0500 

SPATIAL  RESOLUTION  FACTOR  0.0100 

TINE  INTERVAL  FOR  I.B.  TABLE  STORAGE(SEC)  0.1000E-03 

TINE  INTERVAL  FOR  PRESSURE  TABLE  STORAGE  (SEC)  0.1000E-03 

FILE  COUNTERS 

NUMBER  OF  STATIONS  TO  SPECIFY  TUBE  RADIUS  6 

NUMBER  OF  TINES  TO  SPECIFY  PRIMER  DISCHARGE  9 

NUMBER  OF  POSITIONS  TO  SPECIFY  PRIMER  DISCHARGE  6 

NUMBER  OF  ENTRIES  IN  BORE  RESISTANCE  TABLE  7 

NUMBER  OF  ENTRIES  IN  WALL  TEMPERATURE  TABLE  0 

NUMBER  OF  ENTRIES  IN  FORWARD  FILLER  ELEIENT  TABLE  0 

NUMBER  OF  TYPES  OF  PROPELLANTS  1 
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Mini  OF  MM  MTE  DATA  SETS 


2 

Mini  OF  ENTRIES  IN  VOID  FRACTION  TARLECS)  0  0  0 

NUNRER  OF  ENTRIES  IN  PRESSURE  HISTORY  TARLES  8 

NUMER  OF  ENTRIES  IN  REAR  FILLER  ELENENT  TABLE  0 


6ENERAL  PROPERTIES  OF  INITIAL  AMIENT  6AS 


INITIAL  TENPERA^URE  (DEGJl)  529.0 
INITIAL  PRESSURE  (PSI)  U.7 
MLECULAR  HEIGHT  (LBN/LBHOL)  28.896 
RATIO  OF  SPECIFIC  HEATS  1.4000 


GENERAL  PROPERTIES  OF  PROPELLANT  RED 


INITIAL  TENPERATURE  (OEG.R)  529.0 

NININUN  IMPACT  VELOCITY  FOR  EXPLICIT  COMPACTION  HAVE 

(IN/SEC)  10. 

PROPERTIES  OF  PROPELLANT  1 

PROPELLANT  TYPE  JA2  7PF  LOT  472-138 

NASS  OF  PROPELLANT  (LD  17.9000 

DENSITY  OF  PROPELLANT  (LBN/IN**3)  0.0571 

FORM  FUNCTION  INDICATOR  7 

OUTSIDE  DIAMETER  (INS)  0.4150 

INSIDE  DIAMETER  (INS)  0.0400 

LENGTH  (INS)  0.A430 

NUMBER  OF  PERFORATIONS  7. 

SLOT  WIDTH  (HFORI^II)  OR  SCROLL  DIA.  (NF0RN=13)  (INS)  0.0000 
PROPELLANT  STACKED  (0=N0.1=TES)  0 

ATTACHMENT  COIBITION  (0=FREE,1=ATTACHED  TO  TUBE, 

2=ATTACHED  TO  PROJECTILE)  0 

BOND  STRENGTH  (LBF)  (N.B.  ZERO  DEFAULTS  TO  INFINITY)  0.000000 
GRAIN  INHIBITED  ON  OUTER  SURFACE  (0=N0,1=TES)  0 


RHEOLOGICAL  PROPERTIES 


SPEED  OF  COMPRESSION  HAVE  IN  SETTLED  BED  (IN/SEC)  10000. 
SETTLING  POROSITY  1.0000 
SPEED  OF  EXPANSION  IMVE  (IN/SEC)  50000. 
POISSON  RATIO  (-)  0.5000 
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GRAIN  FRACTURE  DATA  FOR  PROPELLANT 

TYPE  1 

NAX.STRESS(PSI) 

BURN  RATE  NULTIPLIER(-) 

DRAG  AND  HEAT  XFER  NULTIPLIER*(-: 

0.000000 

1.00000 

0.000000 

uso.oo 

1.00000 

0.000000 

3625.00 

1.00000 

0.000000 

7250.00 

1.00000 

0.000000 

*  -  ZERO  VALUES 

OF  DRAG  AND  HEAT  XFER  MULTIPLIER'  DEFAULT  TO  BURN  RATE  VALUE 

SOLID  PHASE  THEHMOCHEMISTRT 


NAXINUM  PRESSURE  FOR  BURN  RATE  DATA  (LBF/1N**2)  10000. 
BURNING  RATE  PRE-EXPONENTIAL  FACTOR 

(IN/SEC/PSI*^)  0.4000E-02 
BURNING  RATE  EXPONENT  0.7162 
NAXIMM  PRESSURE  FOR  BURN  RATE  DATA  (LBF/IN**2)  100000. 
BURNING  RATE  PRE-EXPONENTIAL  FACTOR 

(IN/SEC/PSI*^)  0.8881E-0S 
BURNING  RATE  EXPONENT  0.8796 
BURNING  RATE  CONSTANT  (IN/SEC)  0.0000 
IGNITION  TEMPERATURE  (OEG.R)  800.0 
THERMAL  COMXICTIVITY  (LBF/SEC/DEG.R)  0.277QE-01 
THERMAL  OlFFUSIVITY  (1N**2/SEC)  0.1345E-0S 
EMISSIVITY  FACTOR  q.MO 


GAS  PHASE  THERMOCHEMISTRY 


CHEMICAL  ENERGY  RELEASED  IN  BURNIHG(LBF-IN/LBM) 
MXECULAR  UEIGHT  (LBM7LBHOL) 

RATIO  OF  SPECIFIC  HEATS 
CONOLUME 


.20372e«08 

26.8226 

1.2268 

26.9800 


LOCATION  OF  PACKAGEIS) 

PACKAGE  LEFT  BOOY(IHS)  RIGHT  BOOY(INS)  MASS(LBN)  INNER  RAOIUS(IN) 

^  0.000  22.200  17.900  0.000 

PROPERTIES  OF  PRIICR 

CHEMICAL  ENERGY  RELEASED  IN  BURNING(LBF-IN/LHM}  O.OOOOE'KIY 
MOLECULAR  UEIGHT  (LBM/LBMOL)  30.9300 

RATIO  OF  SPECIFIC  NEATS  1.2210 
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OUTER  RADIUS(IH) 
0.000 


SPECIFIC  VDUME  OF  SOLIDdN^ABN) 


23.0000 


PEINER  OISOUME  FUNCTION  (UN/IH/SEC) 


POS.(IHS)  0.00 

TIME(SEC) 

1.00 

3.00 

3.01 

12.00  12.10 

0.000 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.250E-a3 

0.00 

0.00 

0.00 

6.24 

0.00 

0.00 

0.125E-02 

0.00 

0.00 

0.00 

6.24 

6.24 

0.00 

0.150E-02 

0.00 

0.00 

0.00 

0.00 

6.24 

0.00 

0.350E-02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.37SE-02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.500E-02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.525E-02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.725E-02 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

PMUNETERS  TO  SPECIFY  TUBE  QEONETRT 

DISTANCE(IH)  BM>IUS(IN) 


0.000 

2.250 

3.000 

3.090 

19.000 

3.090 

22.000 

2.380 

23.370 

2.360 

208.350 

2.360 

BORE  RESISTANCE 

TABLE 

POSITIONCIHS)  RESISTANCE(PSI) 

0.000 

145. 

0.240 

835. 

0.350 

930. 

0.470 

835. 

0.950 

510. 

3.940 

290. 

187.000 

0. 

THERMAL  PROPERTIES  OF  TUBE 

7.770 
0.2280E-01 
0.700 
529.00 


THERMAL  CONDUCTIVITT  (LBF/SEC/DEG.R} 
THERMAL  DtFFUSIVITT  (IN**2/SEC) 

EMI  SSI  VI TT  FACTOR 
INITIAL  TEMPERATURE  (OEG.R) 
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PROJECTILE  MO  RIFLING  MTR 


INITIAL  POSITION  OF  RASE  OF  PROJECTILEdNS)  22.200 
NASS  OF  PROJECTILE  (LBN)  15.650 
POLAR  NOIENT  OF  INERTIA  (LBO'IN**2)  U.OOO 
ANGLE  OF  RIFLING  (DEG)  0.000 


POSITIONS  FOR  PRESSURE  TABLE  STORAGE 

0.0000  22.2000  30.0000  41.0000  60.0000  90.0000  120.0000  19.0000 

LOCATION  RELATIVE  TO  TUBE  (0)  OR  REAR  OF  AFTERBODY  (1) 

00000000 


TANK  QUN  OPTION  DATA 

NUNBER  OF  DATA  TO  DESCRIBE  AFTERBODY 

TUBE  SURFACE  SOURCE  (0=NONE.1=TABULAR.2=MODELB» 

CENTERLINE  SOURCE  (0=N0NE.l3TABULAR.2=«n0ELED) 

AFTERBODY  SURFACE  SOURCE  (0=NaNE,1-TABULAR,2iWDELE0) 

NUNBER  OF  ENDUALL  DATA  SETS 

NUNBER  OF  PERNEABILITY  DATA  SETS 

NUNBER  OF  REACTIVITY  DATA  SETS 

NUNBER  OF  SEGNENTS  ON  TUBE  SURFACE  SOURCE 

NUNBER  OF  SEGNENTS  ON  CENTERLINE  SURFACE  SOURCE 

NUNBER  OF  SEGNENTS  ON  AFTERBODY  SURFACE  SOURCE 

CONTROL  CHARGE  PRESENT  (O^NO.-I’^YES) 

NUNBER  OF  DATA  TO  DEFINE  EXTERNAL  GEONETRY  OF 
CONTROL  CHARGE  CHANBER 
CONTROL  CHARGE  DETERRED  (0=N0;1=YES) 

BASE  OF  AFTERBODY  NQDELED  AS  DISCONTINUITY 
(0=NO;1=YES) 


GEONETRY  OF  AFTERBODY 

AXIAL  POS.(IN)  RADIUS(IN) 

0.000  0.530 

8.000  1.280 

DESCRIPTION  OF  TUBE  SURFACE  SOURCE 

NUNBER  OF  DATA  TO  DESCRIBE  THICKNESS  OF  LATER 
NUNBER  OF  DATA  TO  DESCRIBE  DENSITY  OF  SEGNEHT  1 
NUNBER  OF  DATA  TO  DESCRIBE  DENSITY  OF  SEGNEHT  2 
NUNBER  OF  DATA  TO  DESCRIBE  DENSITY  OF  SEGHENT  3 


3 

4 
4 
0 
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THIOOKSS  OF  REACTIVE  lAYER 


AXIAL  POS.(IN) 

THICKNESS! IN) 

SECMENT 

0.000 

0.U5 

1 

18.000 

0.145 

2 

22.200 

0.005 

0 

PROPERTIES  OF  SEGMEMT  NUMBER  1 


DEHSITT  OF  REACTIVE  UTER 


DENSITY(UM/IN**3} 

PRESSURE(PSI) 

0.0295 

15. 

0.0455 

11000. 

0.0497 

25000. 

0.0548 

100000. 

NUMBER  OF  DATA  TO  DESCRIBE  BURN 

RATE 

2 

BURN  RATE  DATA 

•MX.PRESS(PSI)  COEFFdN/SEC/PSI*^) 

EXPONENT 

10000. 

0.16000E-03 

1.3010 

100000. 

13.550 

0.0690 

BURN  RATE  ADDITIVE  CONSTANT  (IN/SEC) 

0.0000 

IGNITION  TEMPERATURE  (R) 

800.0 

THERMAL  OOMDUCTIVITT  (LBF/SEC/R) 

0.2770E-81 

THERMAL  DIFFUSIVITT  (IN^2/SEC) 

0.1345E-03 

EMISSIVm  FACTOR(-) 

0.000 

CHBIICAL  EHERGT  (LBF-IN/LBM) 

.93000E+07 

NOLEOXAR  HEIGHT  (LBM/LBMOL) 

22.3900 

RATIO  OF  SPECIFIC  HEATS  (-) 

1.2580 

PROPERTIES  OF  SEGMEHT 

NUMBER  2 

DEIKITT  OF  REACTIVE  LATER 

DENSITY  CLBM/IH^3) 

PRESSURE(PSI) 

0.0295 

15. 
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0.M55 

0.M97 

0.0548 


11000. 

zsooo. 

100000. 


NUMBER  OF  MTA  TO  DESCRIBE  BURN  RATE  1 

BURN  RATE  DATA 

NAX.PRESS(PSI)  0QEFF(IN/SEC/P8I*^)  EXPCHENT 

100000.  0.62520E-04  1.3010 


BURN  RATE  ADDITIVE  CONSTANT  (IN/SEC) 

0.0000 

IGNITION  TEMPERATURE  (R) 

80000.0 

THERMAL  CONDUCTIVITY  (LBF/SEC/R) 

0.277QE-01 

THERMAL  DIFFUSIVITY  (IM^2/SEC) 

0.1345E-03 

ENISSIVITY  FACTOR(-) 

0.000 

CHEMICAL  ENERGY  (LBF-IN/LBM) 

.93000E+07 

MOLECULAR  HEIGHT  (LBM/LBMOL) 

22.3900 

RATIO  OF  SPECIFIC  HEATS  (-) 

1.2580 

BORE  RESISTANCE  DATA  HAVE  BEEN  INTERPRETS)  AS  RELATIVE  TO  PROJECTILE  DISPLACEMENT. 

NOVSUB  MESSAGE...  SETTLING  POROSITY  AT  REFERENCE  COMPOSITION  HAS  BEEN  DEFAULTED  TO 
0.42399  TO  AVOID  IHITIAL  BED  COMPACTION  OF  PROPELLANT  TYPE  1 

VOLUME  OF  AFTERBODY  (IN**3)  21.762 


EQUIVALENT  INTBAL  DATA 


PROJECTILE  TRAVEL(IN) 
CHAMBER  V0LUNE(IH^3) 
GUN  NASS(LBM) 

GUN  RES.  FAC. 

ELEV.  ANGLE(OEG) 
PROJECTILE  MASSLLBM) 


PROJECTILE  TRAV.  (IN) 

0.000 

0.240 

0.350 

0.470 

0.950 


186.900 

622.132 

0.100E^21 

0.000 

0.000 

15.650 


RESISTAHCE(PSI) 

145.000 

835.000 

930.000 

835.000 

510.000 
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3.940 

187.000 


290.000 

0.000 


VEL.  TNRESNOLD  FOR  DYN.  RES.(F/S)  27.000 

VEL.  DEPENDEHCE  OH  OMRCE  UEIGHT(F/S/LaM)  0.000 

ESTINRTED  NUZZLE  VELOCITY(F/S)  0.000 

N.B.  USE  VALUE  FRON  SUOMRT  TABLE.  INTBAL  WILL  NOT  ACCEPT  ZERO 
BORE  AREAdH^)  17.497 

AIR  DEHSITY(LBN/FT^)  0.000 

IGNITER  NASS(LBN)  0.0849 
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